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Abstract 

In this paper a two channel paraunitary filter bank is proposed , which is based on linear canonical transform, instead of discrete 
Fourier transform. Input-output relation for such a filter bank are derived in terms of polyphase matrices and modulation matrices. 
^) It is shown that like conventional filter banks, the LCT based paraunitary filter banks need only one filter to be designed and rest of 
fS| the filters can be obtained from it. It is also shown that LCT based paraunitary filter banks can be designed by using conventional 
power-symmetric filter design in Fourier domain. 
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•"t ' I. Introduction 

O A. Motivation 

Filter banks are now well known tool for time-frequency analysis of a signal HI, lEl- Given a signal x{n), a two 
. channel filter bank (as shown in figlB. splits its Fourier transform (FT) spectrum into two by a low-pass and a high-pass 
C<1 'filter {ho{n), hi{n)} respectively. Output of these filters are then down-sampled to produce analysis filter bank output 
^ ' {uoin) , yi{n)} . To reconstruct the signal {yQ{n),yi{n)} are passed from a synthesis filter bank. The synthesis filter 
^ bank upsamples {yQ{n),yi{n)}, and after passing them from synthesis filters {gQ{n), gi{n)}, adds them to form the 
I filter bank output x{n). The filter bank introduces various distortions like alias distortion, magnitude distortion and 
O phase distortion into the signal. It can be shown that by careful design of the analysis and synthesis filters, all these 
^ distortions could be eliminated and x{n) could be made a delayed version of x{n) by some integer K. Such filter banks 
. . are called perfect reconstruction filter banks (PRFB)lin. |[2l . 

. ^ ' Discrete Fourier transform(DFT), and its generalization z-transform (ZT) find extensive use in the filter bank theory. 
^ Many key theorems of filter bank theory are expressed in Fourier domain or z-domain. Terms like low-pass and high- 
pass make sense in the Fourier domain. Design of the filters is carried out by giving their specifications in the Fourier 
domain. 

In another development, the Fourier transform(FT) has been generalized to fractional Fourier Transform (FrFT), by 
looking at FT as a transform that rotates the signal in time-frequency plane by 7r/2 13], [4]. It is found that the FrFT 
is a special case of three parameter family of transforms called Linear Canonical Transform(LCT)|i5j|, |i6j|. Given a 
continuous time signal x{t), its LCT is given by 



A-„.v.<iM = i/4t /°°i(t)e'(*''-l"'+*"''<!t (1) 

where ad — be = 1. 

In order to define FrFT and LCT on a discrete signal, many sampling theorems in FrFT domain and LCT doamin 
have been derived Q, HI, Q, 111, ifTOl . It is shown Q that if a signal is band-limited in the LCT domain, then it can be 
reconstructed by its uniformly sampled discrete samples. Based on uniform sampling, a discrete time LCT (DTLCT) is 
introduced in fT(F|, and sampling rate conversion theorems for this DTLCT are also derived. 

Since sampling rate conversion is used in the filter banks, it is natural to ask if a sub-band decomposition scheme 
could be developed which is based on the DTLCT. In particular, we want to know how a filter-bank can be designed 
that splits the DTLCT spectrum of an input signal into two sub-band signals. 
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B. Contributions of the paper 

In this paper the samphng rate conversion theorems derived in ifTOl are used to develop a sub-band decomposition 
scheme based on DTLCT. A convolution is defined in the DTLCT domain which is suitable for the sub-band decompo- 
sition scheme proposed. Input-output relation of an LCT filter bank in polyphase and modulation (or alias component) 
domains is derived. Condition for a two channel LCT filter bank to be paraunitary (PU) is also derived. It is shown 
further that a two channel PU LCT filter bank can be designed from a power-symmetric filter in ZT domain. 

C. Notations 

Given a signal x{n), its LCT is denoted by X{uj), FT by X{e^^) and ZT by X{z). Sampling rate of x{n) is denoted 
by Tx- Downsampling x{n) by an integer is denoted by x{n) | N, and upsampling by N is denoted by x{n) | N. 
Matrices are denoted by boldface letters. Complex conjugate of X{uj) is denoted by X^{lo). Conjugate transpose of a 
matrix tl{u>) is denoted by H^f (tj). For a scalar H{z), H^{z~^) is denoted by H{z), and for a matrix H(2;), Yi^{z^^) 
is denoted by H(z). 

II. Review of Filter Banks and DTLCT 

A. Review of DTLCT 

Let x{n) be a discrete signal, obtained by uniform sampling of a continuous signal Xc{t) with a sampling period Tx, 
so that x{n) = Xc{nTx). The DTLCT of x{n), with parameters (a, 5, c, d) is defined as [10|. 

,6/0, (2) 



2 ^nX' - 2nsgn(6)a; + 



where sgn(6) is sign of b. In this paper we will assume 6 > to avoid sgn(6) term. This does not result in loss of 
generality, since all the results in this paper can be obtained for 6 < case in the same way they are obtained for 6 > 
case. 

Upsampling a signal by L is defined as inserting L zeros between two samples. This changes the sampling time of 
the upsampled signal. Thus if y(n) = x{n) | L, then Ty = Tx/ L. It can be shown that iflOj 

Y{uj) = X{Luj) (3) 

Similarly downsampling a signal by M is defined as dropping M — 1 samples out of a block of M samples and 
retaining only the first sample. This also changes the sampling time of the downsampled signal. Thus if y{n) = x{n) [ 
M, then Ty = MTx. It can again be shown that [10] 



^ M-l 

^('^) = T7 E 

m=0 



jdb27rm 

7f;2 [^ + mTT) 

y 



(4) 



B. Review of Filter-banks 

Consider a two channel filter-bank shown in figUl The input-output relation of this filter-bank can be written either 
in terms polyphase matrices, or in terms of modulation (or alias component) matrices | IJ, ||2j. 

Given a discrete signal x{n), having ZT X{z), its type-1 polyphase representation is obtained by writing X{z) as 
X{z) = Xq{z'^) + z-'^Xi{z^), where Xo{z) = EnG/ 3;(2n)z-", Xi{z) = Ene/ ^(2"- + Similarly a type- 

2 polyphase representation is obtained by writing X{z) = Xo{z'^) + zXi{z'^), where Xq{z) is same as above and 
^i(^) = Yl,n&i x{2n — l)z^"-. Using these polyphase representations, it can be shown that 

Xp(z) = Gp(z)Hp(z)Xp(z) (5) 

where 

^^^'^ - X,{z) 
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Gp{z) 



def 



def 



def 



Xi{z) 

Gm{z) 

Hoo{z) 
Hioiz) 



Gioiz) 
Giiiz) 

Hoi{z) 
Hniz) 



(6) 



where the polyphase decompositions of analysis and synthesis filters are of different type (i.e., if {Hq{z), Hi{z)} 
have type-1 polyphase decomposition, then {Go{z),Gi{z)} will have type-2 representation and vice versa), and the 
polyphase decomposition of X{z) and X{z) are of same type. 

The input-output relation can also be expressed in terms of modulation matrices as 



^ni{z) — —Gm{z)'Hm{z)'^m{z) 



(V) 



where 



^m{z) 
Gmiz) 
timiz) 



def 



def 



def 



def 



X{z) 
X{-z) 

X{z) 
X{-z) 

Go{z) 
Goi-z) 

Ho{z) 
Hi{z) 



Giiz) 
Gii-z) 



Hoi 



It can be seen that the polyphase matrices and the modulation matrices are related as 



Hoiz) 
Hi{z) 



Hoi 
Hii 



-z) 



Hooiz) 
Hioiz) 



Hoiiz) 
Huiz) 



(8) 



(9) 



where A; = 1 for type-1 polyphase decomposition and A; = —1 for type-2 polyphase representation. 

A square polynomial matrix H(2;) is said to be paraunitary (PU) if H(2;)H(2;) = dl for some d > 0. It can be seen 
from Q that if polyphase matrix of analysis filters is PU, then its alias component matrix will also be PU. 



If 'tlrn{z) is PU, then by making Gmiz) 



synthesis filters causal, X{z) 



'iHmiz), where K is an integer appropriate enough to make the 



X{z) can be achieved, i.e. the output is delayed version of the input, and perfection 



reconstruction can be achieved. 

If Hm(z) is PU, then it can be seen that Ho{z) satisfies what is called power-symmetry property, i.e. 



Further Hi{z) is related to Hoiz) as 



If Hm(^) is PU, then making 



|F(e^''^)|2 + I//(e-''(^+^))|2 = (i 

Hiiz) = cz'^Hoiz), \c\ = 1, L = odd 
iz) = z^^Hmiz) will give us PRFB. This gives 



(10) 



(11) 



Go{z) = z-^'Hoiz) 

Gi{z) = z-^'Hiiz) (12) 

Thus all the filters of a two channel PU filter bank can be obtained by designing a power symmetric filter Hoiz). 
If Hoiz) is power-symmetric, then Hoiz)Hoiz) is an half-band filter. Thus a power-symmetric filter can be designed 
by designing a half-band filter and computing an appropriate spectral factor (Please refer ||T] for further details.). 
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x(n) 



ho(n) 



Uo(n)^/- ^ yo(n) 



hi(n) 



Ui(n) v_Vyj(n) v_^vi(n) 
Fig. 1 

The two channel filter bank. 



go(n) 



Xo(n) 
x(n) 



g,(n) 



Xj(n) 



III. Convolution and Delay in the DTLCT domain 

Filtering operation in sub-band decomposition is defined as taking product of the DFT(or ZT) of the signal with the 
filter. In time domain this corresponds to convolution between the signal and the filter. In continuous FrFT and LCT 
domain also various convolution theorems have been derived |[Tn . llT2l . 

For LCT sub-band decomposition, it is convenient to define convolution of x{n) and h{n) as 

OO , 2 

/i(n) * a;(n) =^ ^ /i(A;)a;(n - A;)e-^^("-^) (13) 

k=—oo 

Thus if y{n) = h{n) * x{n), then 

j dbm ^ 

Y{uj) = H{uj)X{uj)e'^^ (14) 
With this definition of convolution, it is convenient to define a new delay operator D[.] on a signal x{n) as 

Z)[x](n) =2;(n-l)e'^(-2n+i)^ ^-^g^ 

The properties of D[.] are summarized in the following lemma, 
Lemma 1: The operator D[.] has the following properties, 

1. D^[D''[x]]{n) =D'+^[x](n)foralU,/ G /. 

2. If y(n) = D''[x]{n) then Y{u;) = e-^^'^X{uj). 

3. lfy{n) = x{n) * h{n) then Z)'[x](n) * D^[h]{n) = D^+^[y]{n). 

Proof: The first two results are straight-forward. To prove the last result let 



D\x\{n) * D''[h]{n) = Ll'[x](r)Z)'=[/i]( 



n — r)e b 



r[n—r) 



rel 

= ^x(r-Z)e —b -h{n-r-k)e' '-^-h^^i^-^) (16) 

By re-arranging the index and the power of e, we can obtain D^[x]{n) * D^[h]{n) = D^+\y]{n). ■ 

IV. A Two Channel LCT Filter Bank 

Consider the filter bank in figUl with convolution as defined in (fT3] ). We are interested in writing the input-output 
relation of this filter bank in terms of polyphase and modulation matrices, if ZT (or DFT) is replaced by DTLCT. 

A. Polyphase Domain Analysis 

To define a type-1 polyphase decomposition of x(n), let x{n) = Xe{n) + Xo(ra), where 

, , I x(n) n = even 
xJn) = <^ ^ , , (17) 

^ ' yQ n = odd ^ ^ 

and Xo{n) is similarly defined. Let Xo{n) = D^'^[xo]{n), then X{uj) = Xe{uj) + e^^Xo{uj). Let xo(n) = Xe(n) J, 2 
and xi{n) = Xe{n) [ 2, then type-1 polyphase decomposition of x(n) is given as X{uj) = Xq{2uj) + e^'^ Xi{2oj). 
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x(n)— h (n) -(i 2) -(1 2) go(n) — 



Fig. 2 

A FILTER BANK CHANNEL. 



x(n)^ 



h (n) -(I 2 -y (n) 

(a) 



y (n) 



(b) 



(n) 



- X (n) 



Fie. 3 



Filter bank channels: (a) Analysis filter bank channel (b) Synthesis filter bank channel 



Similarly, we can obtain type-2 polyphase decomposition of x{n) by putting Xoin) = D\xc\{n). In this case X{uj) = 
Xq(2uj) + e~^^ X\(2ijj). Note that the sampling time of the polyphase components T^^ = T^^ = 2Tx. 

The input-output relation of an analysis filter bank channel, shown in fig 3(a) is given in terms of polyphase compo- 
nents as 

Lemma 2: If y{n) = {h{n) * x{n)) [ 2, then 



-jdbuj 



Y{u) 



Xiiu;) 



(18) 



where the polyphase components of h{n) and x{n) are of different type. 
Proof: Let u{n) = h{n) * x{n), then. 



Ue{n) 

Taking DTLCT of both the sides. 



X(.{n) * h(.{n) + Xoin) * ho{n) 

Xe{n) * he{n) + D[xo]{n) * D'^[ho]{n) 

Xe{n) * he{n) + Xo{n) * ho]{n) 



(19) 



- j dhu} 



' {He{uj)Xe{u;) + Ho{u;)Xo{u;)) 

' {Hq{2uj)Xq{2uj) + Hi{2lo)Xi{2uj)) 



(20) 



By noting that Y{u>) = Ue{oj/2), and Ty = 2T„, and substituting it in (l20b . we obtain the required result. ■ 
Similarly the input-output relation of a synthesis filter bank channel, shown in fig |3(b)[ is given in terms of polyphase 
components as- 

Lemma 3: If x{n) = g{n) * (y(n) | 2), then 



" Xoi^) ' 


— jdbuj'^ 

= e '''^y 


' Go(w) ■ 


_ Xi{^) _ 







(21) 



where the polyphase components of both g{n) and x{n) are of same type. 

Proof: Let v{n) = y{n) | 2. Then Xe(n) = Qein) * v{n) and Xo{n) = Qoin) * v{n). This gives D[xo]{n) 
D[go\n * v{n). Taking LCT of both the sides, we have 



Xe{i^) 



-jdbuj 



Go{oj) 



(22) 



By noting that Ty = 2T^, V{io) = Y{2lo), Xe{uj) = Xo(2u;) and so on, and substituting them in ((22]) we get the desired 
result. ■ 
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Using the above two lemma, The filter bank equation for the filter-bank in figHJis 

Xp(u;) = e 



- j dhu> 



Gp(w)Hp(a;)Xp(u;) 



(23) 



where 



Xp(a;) 



Xp(a;) 



H(a;) 



Xo(a;) 






Xi{u^) _ 




















Gio 






Gil 








(a;) 




^11 


(^) 



(24) 



5. Modulation Domain Analysis 



To write the filter bank equation in terms of alias component matrices, it is preferable to consider a channel of the 
filter bank as shown in fig|2]. The input-output relation for this channel is given as foUows- 
Lemma 4: For the channel in figl2], the input and output are related by 



Moreover, X{lu + tt) is given by 

X{L0 + 7r) 



1 



1 —jdb{uj-\-m7r)^ 

m=0 



G{uj)H{lo + m'Tr)X{uj + mvr) 



1 —jdb{LJ-\-m7r)^ 

- ^ e ^ G{uj + 7t)H{u + rmT)X{uj + mvr) 



m=0 



Proof: By using the following relations 



[/(a;) 

Viu;) 
X(uj) 



' F{uj)X{uj) 



1 



1 jdb2rmT{uJ + rmr) uj 



T2 

m=0 2/ 

y(2a;) T^=Ty/2 = T^ 



jdhuj 
2T'i 



G{u)V{uj) 



and some algebraic manipulation we can obtain (125]) . 
To prove (l26l) . we note that using (|25] ). we have 



^(^ + vr) = -E 
Further noting that for any X{uj) 



1 -jdi)(^ + (m + l)7r)^ 



'G{uj + 7r)i?(w + (m + l)7r)X(a; + (m + l)7r) 



J dEi27r{ijj + 7r) 

+ 27r) = e 5^ 



and substituting from (|29ll into ([28]) we get ( |26l ). 

Using the above lemma we can write the following input-output relation for the filter bank shown in figH] . 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 
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Xm(a;) 
Gm(a;) 
Hm(w) 



X(w + 7r) 

l(u; + 7r) ^ 

Go (a;) Gi{uj) 
Go(a; + 7r) Gi(a; + 7r) 



e 








_ jdb(c£+7ryf 



(31) 



Using the relation X{uj + 2tt) 



where 



■ 2nd b I 



X{ijj), It can be shown that 

X„(c^) = B(u;)C(e^'")Xp(u;) 



(32) 



B(^) 

CM 



1 

6 ^2 

1 -eJ"^ 



i(2a;+7r) 



(33) 



V. Two Channel Paraunitary Filter Banks 



A square matrix 'il{oj) is defined to be a paraunitary matrix if HJ^(w)H(u;) = dl for some d > and for all Let 
Hm(tj) in (I30ll be PU, then by setting 

Gm{oj) = e-^'^" A, (a;)H^, (l^) (34) 

where K is an integer (needed to make the synthesis filters causal), we can obtain perfect reconstruction. 

We now show that as in the case of ZT, H(u;) can be made PU if Hq{uj) is power-symmetric, and Hi{uj) can be 
obtained from Ho{lo) by inspection. 

Lemma 5: If Hq{uj) is power symmetric, i.e. 



\Ho{uj)f + \Ho{uj + Tr)f 



1, 



and Hiiuj) is such that 



H 



^[uj)=e^y ^ ' I Hq^[uj + -K), 



(35) 



(36) 



for an integer L, then H(u;) is PU. 

Proof: By expanding each element in Hm(w)H^^(t(;),substituting from (1351) and (l36l ). and using the relation 

- 2'Kdb / \ \ m 

H,{uj + 27r) = e^~?^^'^^''' Hi{oj), where f G {0, 1}, we can get H„(a;)H^,(a;) = /. ■ 
The above lemma suggests that just like in the ZT filter bank, problem of designing a two channel PU LCT filter bank 
can be reduced to the problem of designing a power-symmetric filter, and rest of the filters could be obtained from this 
filter. However a power-symmetric filter in LCT domain can be obtained from a power-symmetric filter in ZT domain, 
as shown in the following lemma. 

Lemma 6: Let Hq{uj) = DTLCT{ho). Let ho{n) be defined as ho{n) = ho{n)e^i^ . Let Ho{e^'^) = FT{ho). If 
Ho{e^^) is power-symmetric, then Hq{lo) will also be power-symmetric. 
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Fig. 4 

Magnitude of the input signal X{uj). 



Proof: Hq{uj) can be written as 



jdbuj 



From this it can be seen that 



(37) 



Ho{u;)H*,{u) = H^{enm{e-n (38) 
and so if HQ{e^'^) is power-symmetric, then Hq{lo) will also be power-symmetric. ■ 

VI. An Example 

In this section we take an example which serves as a demonstration of the LCT filter bank. In this example we take 
the LCT to be FrFT at an angle tt/A. The [a, b, c, d] matrix will thus be 



a b 




c d 





cos(7r/4) sin(7r/4) 
— sin(7r/4) cos(7r/4) 



(39) 



The input signal x{n) is such that its LCT X{lo) has peaks at{|2|, ^} , as shown in figH The samphng 

time of the input signal x{n) is taken as T = 0.05. 

An LCT power-symmetric filter is generated by applying lemma [6] on a ZT power-symmetric filter. In this example, 
the ZT power symmetric filter h{n) is taken from LI l,example 5.3.2. The magnitude response is shown in fig©. Using 
lemma m the LCT power symmetric filter ho{n) is given by 



/io(n) = h{n)e 2 b 
Filter hi{n) is obtained from hQ{n) by applying lemma|5] This gives /ii(n) as 

hi{k) = ho,{N - k)e^'f'^''\k e {0..N} 

where 



{N - kf + k^ +{N - k)TT 



db 
2r2 



TT 



(40) 

(41) 
(42) 
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Fig. 5 

MAGNITUDE(dB) of the power symmetric filter H{e^'^), THAT IS USED TO GENERATE Ho{u}) 




-3-2-10123 -3-2-10123 



(a) (b) 

Fig. 6 

MAGNITUDE(DB) of the LCT of various filters USED IN THE LCT FILTER BANK: (A) ANALYSIS FILTERS 

{i?o(w), -ffi(w)} (B) Synthesis filters {Go{u}),Gi{ui)} 



The magnitude response of the LCT of analysis filters is shown in fig |6(a)| 

Finally the synthesis filters {Go{io), Gi{uj)} are obtained from {Ho{lu), Hi^lo)} using (l34l ). which gives 

g,{k) = h^N - A:)e-2^[(^-'=)'+'='U G {O..N},i G 0,1 (43) 
The magnitude response of the LCT of synthesis filters is shown in fig |6(b)| 

The LCT of analysis filter output {yQ{n),yi(n)} are shown in fig |7(a)| and fig |7(b)| The peaks at {^^f , ^5^} in 
Yo{uj) correspond to peaks at {|^, in X{u>). Similarly peaks at , in Yi{uj) correspond to peaks at 

f 4827r 4127r\ ■ yf, A 

iw -512 f in 

The signal is then reconstructed by passing {yo{n'),yi{y)} from synthesis filter. The LCT of the synthesis filter 
output X{u!) is shown in figlH It is observed that the output matches with input delayed by N samples. 
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(a) 



(b) 



Fig. 7 

Magnitude of analysis filter bank output: (a) Fo(tj) (b) YiIuj) 
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0.3 - 




-3-2-10123 



Omega 

Fig. 8 

Magnitude of the filter bank output Xiuj). 



VII. Conclusion 

In this paper a two channel paraunitary fiher bank based on Unear canonical transform has been developed. This kind 
of fiher bank can be useful in sub-band decomposition of the signals that are not band limited in the Fourier domain, but 
band limited in an LCT domain. Input-output relations for such a filter bank in polyphase and modulation domain are 
derived. It is also shown that such a filter bank can be designed by using standard design procedure for power-symmetric 
filters in Fourier domain. 

The future work in this direction may include generalization of LCT based filter bank to more than two channels, and 
development of LCT based cosine modulated filter banks. 
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